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Abstract: We have previously demonstrated a pulsed direct detection IPDA lidar to 
measure range and the column concentration of atmospheric CO 2 . The lidar measures the 
atmospheric backscatter profiles and samples the shape of the 1,572.33 nm CO 2 absorption 
line. We participated in the ASCENDS science flights on the NASA DC-8 aircraft during 
August 2011 and report here lidar measurements made on four flights over a variety of 
surface and cloud conditions near the US. These included over a stratus cloud deck over 
the Pacific Ocean, to a dry lake bed surrounded by mountains in Nevada, to a desert area 
with a coal-fired power plant, and from the Rocky Mountains to Iowa, with segments with 
both cumulus and cirrus clouds. Most flights were to altitudes >12 km and had 5-6 altitude 
steps. Analyses show the retrievals of lidar range, CO 2 column absorption, and CO 2 mixing 
ratio worked well when measuring over topography with rapidly changing height and 
reflectivity, through thin clouds, between cumulus clouds, and to stratus cloud tops. The 
retrievals shows the decrease in column CO 2 due to growing vegetation when flying over 
Iowa cropland as well as a sudden increase in CO 2 concentration near a coal-fired power 
plant. For regions where the CO 2 concentration was relatively constant, the measured CO 2 
absorption lineshape (averaged for 50 s) matched the predicted shapes to better than 1% 
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RMS error. For 10 s averaging, the scatter in the retrievals was typically 2-3 ppm and was 
limited by the received signal photon count. Retrievals were made using atmospheric 
parameters from both an atmospheric model and from in situ temperature and pressure 
from the aircraft. The retrievals had no free parameters and did not use empirical 
adjustments, and >70% of the measurements passed screening and were used in analysis. 
The differences between the lidar-measured retrievals and in situ measured average CO 2 
column concentrations were <1.4 ppm for flight measurement altitudes >6 km. 

Key words: atmospheric CO 2 ; IPDA lidar; CO 2 DIAL lidar 


1. Introduction 

Although increasing atmospheric CO 2 is widely accepted as the largest anthropogenic factor 
causing climate change, there is considerable uncertainty about its global budget. To be useful in 
reducing uncertainties about carbon sources and sinks, studies have reported that space-based 
atmospheric CO 2 measurements need -0.3% precision and stability, on regional scales of from 10 
degrees square [1-3] to 1-5 degrees square [4], Several groups have analyzed space missions using 
passive spectrometers [5-8], and the GOSAT mission [9] is now making global CO 2 measurements 
from space using an FTIR spectrometer and reflected sunlight. However, limitations from sun angles, 
atmospheric scattering, and detector dynamic range restrict its measurements to daytime cloud-free 
scenes in the low and mid-latitudes. An inherent error source with space-based passive spectrometers 
is optical scattering from aerosols and thin clouds, particularly cirrus, in the illumination or 
observation paths [10,11], Even small amounts of scattering in either of these paths modifies the 
optical path length and thus the total CO 2 absorption measured, and so it can cause large retrieval errors 
even for thin cirrus clouds [11], CO 2 absorption measurements made at slant angles using km-sized 
footprints are also susceptible to errors caused by changes in the optical path length caused by 
variability of the surface reflectivity within the footprint. 

To overcome these limitations, the US National Research Council’s 2007 Decadal Survey 
recommended a new space-based CO 2 measuring mission called ASCENDS [12] using the laser 
absorption spectroscopy approach. The mission’s goals are to quantify global spatial distribution of 
atmospheric CO 2 with ~1 ppm accuracy, and quantify the global spatial distribution of terrestrial and 
oceanic sources and sinks of CO 2 with approximately monthly time resolution. The ASCENDS 
approach directly measures range and offers column CO 2 measurements made through thin clouds and 
aerosols, continuous measurements over the cloud-free oceans, at low sun angles and in darkness, all 
of which are major advantages over passive sensors. The ASCENDS mission organizers held a 
workshop in 2008 to better define the science and measurement needs and plans for future work [4], 
ESA has also conducted mission definition studies for a similar space mission called A-SCOPE [3,13], 
and their lidar sensitivity and spectroscopic analyses have been published [14,15], 


Remote Sens. 2014, 6 


445 


2. CO 2 Line Choice and Lidar Approach 

The IPDA lidar technique is widely used for open-path laser absorption spectroscopy 
measurements [16,17]. Typically two laser wavelengths are used, and the target is illuminated with the 
laser alternatively tuned onto and off the gas absorption line. The path-integrated gas absorption 
attenuates the on-line laser energy relative to the off-line wavelength. If one knows the difference in 
gas absorption cross-sections at the on- and off-line wavelengths, the elevations of the measurement 
endpoints, and the atmospheric temperature and pressure in the path, and measures the range and the 
optical depth of the gas absorption line, one can solve for the path-integrated gas number density and 
the column concentration. Several groups have made airborne measurements of atmospheric CO 2 
using IPDA lidar, and summaries have been recently published [18-20], 

Our airborne lidar [21-24] was developed to demonstrate a dual-channel pulsed IPDA approach as 
a candidate for the ASCENDS mission, and a description of its configuration and its performance 
during flights in 2009 has been recently summarized [25], The approach uses two tunable pulsed laser 
transmitters allowing simultaneous measurement of the absorption from a CO 2 absorption line in the 
1,570 nm band, O 2 extinction in the oxygen A-band, and surface height and atmospheric backscatter in 
the same path. A tunable laser is stepped in wavelength across a single CO 2 line for the CO 2 column 
measurement, while simultaneously a laser is stepped across a line doublet near 765 nm in the Oxygen 
A-band for an atmospheric pressure measurement [26-28]. Both lasers are pulsed at a 10 kHz rate, and 
the two absorption line regions are repeatedly sampled in wavelength steps at -300 Hz. Both laser 
trans m itters utilize tunable diode lasers followed by laser amplifiers. The direct detection receivers 
measure the time resolved laser backscatter from the atmosphere along with the energies of the laser 
echoes from the surface. After suitable averaging the gas extinction and column densities for the CO 2 
and O 2 gases are estimated from the sampled wavelengths of the surface reflected line shapes. 

We use a single line in the 1,570 nm band for the CO 2 measurement [10]. This vibration-rotation 
band of CO 2 has an appropriate range of absorption that provides good sensitivity to the surface echo 
signal and to variation in CO 2 in the lower troposphere. This band has minimal interference from other 
atmospheric species like H 2 O, and has several temperature insensitive lines. The shorter wavelength 
lines in the R-branch are a better match to available laser and detector technologies. The R16 line at 
1,572.335 nm has been analyzed and is attractive for CO 2 measurements [29]. It has low temperature 
sensitivity, particularly to changes in the lower atmosphere. 

Our approach samples the CO 2 line shape at multiple wavelengths and we used 30 samples for 
these flights. Sampling at multiple wavelengths across the absorption line provides several capabilities 
including the use of using several atmospheric weighting functions in retrievals [29], This approach 
also provides information that allows solving for several different instrument parameters, such as 
baseline tilts, wavelength offsets, and wavelength dependent lidar transmissions. Our work has shown 
this information is quite important to minimize instrumental biases in retrievals. It also allows solving 
for useful spectroscopic information, such as line center wavelengths, line widths and errors in the fits [30]. 

Using pulsed lasers and receiver processing to time (and height) resolve the laser backscatter profiles 
also allows post detection signal processing to isolate the laser echo signals from the primary scattering 
surface, and to reject backscatter from the atmosphere that arrives earlier. Hence it allows isolating the 
full column measurement from potential bias errors caused by atmospheric scattering [10,11], Isolating 
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the surface echo pulse with a time gate in receiver signal processing also substantially improves the 
receiver’s signal-to-noise ratio (SNR) by limiting the amount of noise from the detector and the solar 
background. For space missions, the photon efficiency of the lidar measurement technique is quite 
important, and the photon efficiencies for three IPDA lidar approaches have been recently analyzed and 
compared [31]. 

3. Airborne Lidar Description 

For these flights on the NASA DC-8 aircraft [32] we adapted the previous version of the airborne 
lidar [25] flown on the NASA Glenn Lear-25 aircraft in 2009. A photograph of the flight configuration 
is shown in Figure 1. The sensor head was mounted above the aircraft’s aft-most nadir window 
(Port 9). The window assembly uses separate wedged and anti-reflection coated optical windows for 
both the transmitter and received beams to minimize optical crosstalk and etalon fringing. The airborne 
lidar specifications are listed in Table 1. 

Figure 1. (a) NASA DC-8 aircraft. CO 2 Sounder instrument installed above Port 9 on the 
DC-8, (b) Rack with CO 2 lidar electro-optics, (c) Enclosure with transmit optics and 
receiver telescope that is coupled, via fiber optics, to racks, (d) Rack with O 2 lidar electro- 
optics in the telescope side of the rack, with flight computers on the opposite side for 
operators, and Picarro in situ gas analyzer. 



The approach used to sample the CO 2 line and the lidar’ s block diagram is shown in Figure 2. The 
seed source for the laser signal is a single frequency DFB laser diode, which is operated near 1,572.33 
nm by controlling its temperature and current. A ramp from a signal generator was used to sweep the 
current to the diode laser, and hence its output wavelength. The diode’s CW output is then gated into 
pulses using an acousto-optic modulator (AOM) to an Erbium Doped Fiber Amplifier (EDFA). For 
these measurements, 30 wavelength samples are used across the CO 2 line. A small percentage of the 
CW seed laser output is split off and directed through a fiber-coupled CO 2 absorption cell and to a PIN 
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detector. The C0 2 cell serves as a monitor for center wavelength of the sweep. The laser output is a 
sequence of 1 us-wide laser pulses that occur at a 10 KHz rate and the collimated transmitted laser 
signal exits towards nadir though the transmit window. 


Table 1. 201 1 Pulsed Airborne C0 2 Lidar Parameters. 


C02 line 

C02 line center wavelength 
Laser min & max wavelengths 
Laser wavelength steps across line 
Laser wavelength change/step 
Laser peak power, pulse width, energy 
Laser divergence angle 
Seed laser diode type 
Wavelength monitor 
Calibration heterodyne detector 
Laser Pulse Modulator (AOM) 

Fiber coupled C0 2 cell 
Fiber Laser Amplifier (EDFA) 

Laser line scan rate 
Laser linewidth for each step 
Receiver Telescope type 
Telescope diameter 
Receiver FOV diameter 
Receiver optical bandwidth 
Receiver Optics Transmission 
Detector PMT type 
Detector quantum efficiency 
Detector dark count rate 
Receiver signal processing 
Histogram time bin width 
Receiver integration time 
Recording duty cycle 


R16, 6,359.96 cm -1 
1,572.335 nm 

1,572.28 nm, 1,572.390 nm 
30 

~ 3.8 pm (0.0154 cm 4 ) 

25 watts, 1 jlis, 25 jliJ 
100 jurad 

DFB: Fitel FOL15DCWD 
Burleigh WA1650 
New Focus 2053 -FC InGaAs 
NEOS Model: 26035-2-155 
80 cm path, ~200 Torr pressure 
IPG EAR- 1 OK- 1571 -LP-SF 
300 Hz 
~15 MHz 

Cassegrain, f/10 (Vixen) 

20 cm 
200 prad 
800 pm FWHM 
0.64 

Hamamatsu H10330A-75 
4% (these flights) 

~ 500 kHz 

Photon counting/histogramming 
8 ns 

0.9 s per readout 
90% (0.9 s every 1 s) 


Figure 2. (Left) Illustration of wavelength sampling approach, that samples the 1,572.33 nm 
C0 2 absorption line at 30 wavelengths at a 300 Hz rate. The lidar parameters are 
summarized in Table 1. (Right) Block diagram of the C0 2 channel of the airborne lidar. 
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The output laser wavelengths during the scan were precisely calibrated using a laser heterodyne 
technique several times during most flights in the campaign. The technique is described in Appendix A. 
The analysis showed that modeling the wavelength scan using 5 contiguous piecewise quadratic 
functions, with the slopes matched at their endpoints, would reproduce the scan wavelengths to within 
0.1 pm. For the data analysis for each flight, the fits were made to calibration samples taken during 
that flight. If a calibration was not performed during a flight (e.g., for Iowa), data from the previous 
wavelength calibration was used. Analysis of the scan shape taken on different days showed little 
change in wavelength spacing, particularly around the CO 2 absorption line peak. 

The laser backscatter is collected by the receiver’s 20 cm diameter Cassegrain telescope, which 
views nadir through the receiver window. A multimode fiber is used to couple the optical signal from 
the telescope focal plane to the receiver optics. After passing through an optical bandpass filter, the 
signal is focused onto a PMT detector. The electrical pulse output from the PMT was amplified and 
passed through a threshold detector. The gain of the PMT used in the 2011 flights allowed a single 
photon counting efficiency of ~4% near 1,572 nm. Typical PMT dark count rates were 500 KHz and 
when viewing sunlit desert surfaces the total noise count rate from the detector was a few MHz. The 
pulses from the receiver’s PMT detector and discriminator are binned and accumulated by the 
multi-channel scalar (MCS). The MCS had a time resolution of 8 ns/bin that allowed ~5 receiver 
samples on both the leading and trailing edges of the laser pulses. The start time of the MCS sweep is 
synchronized with the trigger for the first laser pulse of the 30, and hence the start of the pulsed 
wavelength sweep. Each MCS sweep contains a time-resolved histogram of PMT pulse counts vs. 
time. This contains the complete time resolved lidar backscatter record from the laser firing to the 
ground for all 30 pulses. At the end of 0.9 s, each MCS bin contains the total receiver counts for the 
300 laser sweeps. Although the signal strength varied with surface reflectivity and altitude, typical 
received signals for off-line wavelengths were -9000 counts for each laser pulse. The receiver 
histogram record is then read and stored. The laser trigger and data acquisition are synchronized to 
timing markers from the GPS receiver and data was stored every second. The computer also digitizes 
other signals, including the inertial guidance system output from the aircraft and GPS position and 
time. A nadir viewing video camera also records the visible image though the nadir window. 

4. CO 2 Measurement Processing and Retrievals 

For these flights the lidar recorded the time- and wavelength-resolved laser backscatter for 300 line 
sweeps during each one-second recording interval. After the flight for each record we used 
measurement analysis approach shown in Figure 3 to process each one-second record. The range to the 
reflecting surface is determined using the approach described in Section 7. Figure 3 shows a typical 
record with an elevated cloud at a range of 5.7 km and the ground at 7.6 km. The range and pulse 
energies at all wavelengths can be extracted for either the strongest returns, for those nearest the 
predicted ground level, or for ground along with two cloud layers. Since the receiver records the entire 
backscatter profile at all wavelengths, this also permits extracting the wavelength-averaged range-resolved 
backscatter profile at the off-line wavelengths. 
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Figure 3. Measurement processing approach, using examples from lidar receiver 
measurements with 1 s integration time from the flight over Iowa, (a) The 30 wavelength 
stepped laser pulses reflecting from the surface are evident in the time averaged 
backscatter profile, plotted with 8 ns/MCS bin. (b) A time-expanded view of a sample 
backscatter profile recorded for an offline pulse integrated over 1 s. This sample shows the 
reflected signal from a thin cloud at a distance of 5.7 km from the aircraft and the echo 
pulse from the ground at a range of 7.6 km. (c) Time expanded view of two pulse 
reflections (overlayed) from the surface, with red an offline pulse, and blue a pulse at a 
wavelength near the CO 2 line absorption peak, (d) The photon counts for each pulse are 
computed by summing the counts between the pulse edges, and with the wavelength scan 
calibration give the transmission shape for the CO 2 line. Although the atmospheric 
scattering above the ground causes a loss of energy, the measured shape of the CO 2 line in 
the column to the ground is otherwise unaffected. 
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We followed the extraction step by the CO 2 processing and retrieval algorithm shown in Figure 4. 
As part of the data processing step, the recorded data is tested against criteria and is screened using the 
algorithm described in Appendix B. For all flights, >70% of the measurements passed the screening 
criteria. The algorithm averages the extracted values and solves and corrects for instrumental offsets. 
For these flights these solved for a wavelength offset and for a slope in the lidar’ s wavelength 
response. The algorithm then calculated the expected line shape from an atmospheric profile. This was 
based on either vertically resolved data from the in situ measurements during the aircraft’s spiral down 
maneuver (DC-8 atmosphere), or from vertically resolved data extracted from the MERRA atmospheric 
model at the nearest time and location to those of the measurement (MERRA atmosphere). These 
sources are described more fully in Section 6. The line shape calculation also accounts for the nadir angle 
of the aircraft, so also allows analysis of off-nadir measurements, such as those made in a banked turn. 
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Figure 4. Data analysis and retrieval approach used to fit the CO 2 lineshape and determine 
column average CO 2 abundance and other parameters. The MCS record of the backscatter 
profile is analyzed to identify the pulse edges, solve for range, and determine the integrated 
counts as shown in Figure 3. These are then normalized by the transmit energy monitor 
data to give the measured lineshape. In parallel, the algorithm computes the predicted 
lineshape based on the optical path traveled by the lidar pulses, the atmospheric conditions 
of the time. The vertically resolved atmospheric state is computed using either in situ 
measurements from the aircraft’s spiral down manuever or data from the MERRA 
atmospheric model. Initially the concentration of 390 ppm was used at all altitudes. The 
measured and predicted line shapes were independently averaged over 10 s, the 
measurements and calculations were compared, and the residual was computed. The 
algorithm then adjusts parameters to optimize the fit. It first varied the wavelength offset 
and instrument baseline slope (left hand side) to adjust measurements and then it varied 
the CO 2 column concentration (right hand side) in the prediction. The result is the 
retrieved CO 2 column concentration that minimizes the error. The RMS error in the fit, the 
confidence interval of the measurement, and OD values are also retrieved and are used for 
screening out poor fits and bad data, as described in Appendix B. 
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The spectroscopic parameters for the CO 2 line were from HITRAN 2008 [33] and used a line 
mixing correction [34]. The column average CO 2 concentration in the model atmosphere is then 
adjusted to minimize the difference between it and the lidar observations. An example of the line fit 
process is shown in Figure 5, and an error analysis for the lidar measurements is summarized in [25]. 

Using this approach to simultaneously determine CO 2 concentration and instrument model 
parameters considerably improves the measurement accuracy by solving and correcting for slowly 
varying instrument offsets and nonlinearities. After the fit several parameters are extracted including 
path averaged concentration, the measured line shape, the line fit’s confidence limit, and the rms error 
in the fit. The line fits were usually good with rms errors of <1% and typically less than twice the 
signal shot noise limit. From the fitted line we also compute and report the difference in the optical 
depths of the peak wavelength and ±50 pm offset points, that is 

DOD(pk, 50) = OD(A pk )-{OD(A pk - 50 pm)+OD{A pk + 50 pm)} / 2 (i) 

where the OD values were measured on the fitted line shapes. Examples of fitted line shapes measured 
at different altitudes and these points are shown in Figure 6. 


Figure 5. An example line shape comparison before and after the concentration 
adjustment. A sample 10 s averaged transmission CO 2 lineshape measured (blue dots) for 
the Iowa flight from 4.4 km altitude. The initial comparison is to a calculated lineshape 
with 390 ppm concentration that gives a noticeable fit error (blue circles in top plot). After 
optimizing the XCO 2 in the calculations to 378 ppm, there is smaller residual (green circles 
in top plot) and rms error. The difference between the initial calculations (red line in 
bottom plot) and the optimized fit for 378 ppm (green line) is small when compared to the 
overall transmission lineshape. 
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Figure 6. Example of fitted C0 2 line shapes, plotted as optical depths, for previous 
airborne lidar measurements at the altitudes indicated. The optical depths at the fitted line 
peak and the average of the fitted values at the sides, at peak wavelength ±50 pm, are used 
to compute the value of DOD(pk,50) in Equation (1) shown in subsequent plots. 



5. 2011 Airborne Campaigns 

During August 2011 we flew an ASCENDS airborne campaign, and after a series of test and 
checkout flights, four science flights were made. Measurements were made from stepped altitudes 
from 3 to 13 km over locations including a stratus cloud deck over the Pacific Ocean, two desert areas 
in the US southwest, and over the US Midwest from the Rocky Mountains to Iowa. The flight 
locations, flight patterns and atmospheric conditions are summarized in Table 2 and diagrams of the 
flight’s ground track patterns are shown in Figure 7. 


Table 2. Summary of 201 1 Flights reported here. 


Notation: 

Flight 2 

Flight 3 

Flight 5 

Flight 6 

Location: 

Pacific Ocean 

Railroad Valley, NV 

Four Corners, NM 

Iowa 

Flight Dates (all 2011) 

2 August 

3 August 

9 August 

10 August 

DC-8 Spiral Down Location (Latt, Long) 

33°02'N 

38°34'N 

36°44'N 

41°43'N 

122°57'W 

1 15°47'W 

108°13'W 

9 1 °2 1 'W 

Surface Elevation (m) 

0* 

1400 

700 

235 

Surface Pressure 1 (mbar) 

1013 

850 

924 

988 

Approx Duration (min) 

264 

278 

333 

390 

Primary Altitude range (km) 

3.2-12.8 

2.9-11.2 

4.5-12.2 

3-12.1 

Altitude step size 2 (km) 

3 

1.7 

3 

1.5 

Number of altitude steps 3 

3-5 

6-8 

3M 

8-10 

Time span of data analyzed (UTC hours) 

16.25, 18.3 

19.2,21.9 

17.3,20.25 

21.0, 23.25 


Notes: 1 Surface Pressure at spiral location (is listed as ocean surface for Flight 2), * elevation of ocean surface; 2 Altitude 
step sizes varied widely between and within flights; 3 Number of altitude steps: first number indicates the number of steps 
after the spiral to the highest altitude. Second number indicates the total number of steps and includes any steps done 
before the spiral and after the highest altitude following the spiral. 
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These flights allowed testing the performance of the improved lidar and optimizing the retrieval 
algorithms. As part of this, the lidar’ s measurement of range, concentration and differential optical 
depth were made under a wide variety of conditions. Those were compared against predictions made 
using data from atmospheric model and the in situ sensors. 


Figure 7. Ground track patterns for the flights described in this paper, (a) Pacific Ocean 
west of Baja California, (b) Railroad Valley NV and surrounding mountains, (c) near Four 
Corners NM, and (d) near Iowa City Iowa. The flight patterns were based on surface 
topography, air traffic and other considerations. 




Longitude 




Longitude Longitude 


6. Airborne C0 2 Measurements and Calculations 

The NASA DC-8’s REVEAL instrumentation included a GPS receiver and in situ sensors for air 
pressure, temperature and water vapor. The LaRC AVOCET in situ C0 2 sensor [35,36] was also part 
of the campaign and flew on the DC-8. We used readings from AVOCET gas analyzer, which reported 
data at a 1 Hz rate, for all comparisons of column CO 2 concentrations. The altitude distribution of the 
CO 2 concentrations measured by the AVOCET sensor is plotted in Figure 8 for the four flights. The 
plots show the concentrations for most flights and altitudes were near 390 ppm. There are several 
exceptions, including the lower concentration for Flight 6, especially below 2 km. The CO 2 spike for 
the power plant plume for Flight 5 is also notable, as is the lower concentration below 1 km for Flight 2. 
The column-integrated water vapor mixing ratio is also plotted in Figure 8 for all flights. The results 
show that for the air columns measured in the flights the water vapor mixing ratios varied from 1% to 
2% near the surface and in the boundary layers, but generally converged to 0.4% to 0.5% for columns 
measured from aircraft altitudes near 12 km. 

In order to compute the expected line shapes and retrievals based on the lidar measurements, we 
used two different sources to estimate the vertical profiles of atmospheric temperature, moisture and 
pressure beneath the aircraft. One source was the DC-8’s REVEAL system that reported temperature, 
pressure and water vapor in situ readings at a 1 Hz rate. All flights had one spiral down maneuver from 
the highest altitude to near the surface permitting the in situ sensors to sample the entire column almost 
to the surface near the center of the flight pattern. Retrievals based on line shapes computed from the 
DC-8’s REVEAL sensor readings are the most accurate but only were available for the column at the 
spiral down locations. 
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Figure 8. (a) In situ CO 2 concentration vertical profile measurements from AVOCET in 
the spiral-down segments of the four flights in August 2011 over the Pacific Ocean (in 
blue), in Railroad Valley, Nevada (in black), the Four Comers, NM (in red), and West 
Branch, Iowa (in green), (b) The column water vapor mixing ratios as a function of flight 
altitude above sea-level from the DC-8 in situ water vapor measurements for the four 
flights, measured during their the spiral down segment. 


DC-8 AVOCET Doto, 2011 flights 




(b) 


The other approach to compute line shapes sampled values from a layered atmospheric model for 
the nearest location and time of day. For this we used the meteorological analysis data from the 
Goddard Modeling and Assimilation Office, Modem-Era Retrospective Analysis For Research and 
Applications (MERRA) [37], Data were drawn from the 4/day (00, 06, 12, 18 UT) analyzed fields on 
the full model with 72 vertical layers on a grid of 0.5 degree latitude x 0.67 degree longitude. We 
extracted the latitude-longitude interpolated MERRA data at the nearest positions in time along the 
DC-8 flight track. The parameters extracted include surface pressure, and atmospheric pressure, 
geopotential height, temperature, water vapor specific humidity, and winds for the lowest 42 layers, which 
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were sufficient to cover the DC-8 flight altitudes. Since MERRA models the atmosphere for the entire 
globe, this approach also allows a preview of possible retrieval results from a satellite-based lidar mission. 

Figure 9 shows for the computed DOD(pk,50) vs*, vertical distance above the measurement surface, 
with computations based on the in situ readings. To aid their comparison the plots are for the 
DOD(pk,50) values for the flights minus those for a US standard atmosphere with a constant mixing 
ratio of 390 ppm. At the same distance above the surface, the DOD values differ by several percent for 
the different flights. Due to the summer conditions, the air is warmer and the computed DOD(pk,50) 
values for the flights are lower than those of the US standard atmosphere. The colder air over the 
Pacific Ocean had the highest number density and computed DOD values. 

Figure 9. Plot of DOD(pk,50) calculated for the four flights from AVOCET and DC-8 
in situ measurements as a function of range to surface. In order to highlight the differences, 
the same calculation for a US Standard Atmosphere with a vertically uniform CO 2 
concentration of 390 ppm was subtracted. The colder temperatures for flight 2 (over the 
Pacific Ocean) increased its air density and resulted in a lower difference with the US 
Standard atmosphere. 



7. Range Measurements 

Determining the value of the column number density and concentrations from an IPDA lidar 
requires an accurate measurement of range [25]. During the 1 s over which these measurements were 
integrated, the DC-8 aircraft travelled 200-250 m horizontally. Over flat terrain the range is usually 
uniform and the recorded echo pulse, averaged over 1 s, maintained the shape of the transmitted laser 
pulse. However, variations in the elevation of the reflecting surface over this distance (from rougher 
ground or cloud tops) as well as aircraft pointing and attitude changes can vary the range during the 
integration time. These can cause a broadening of the averaged echo pulses. 

For previous work we used a cross correlation algorithm to estimate range [24]. For this analysis we 
used timing to the echo pulse’s centroid, which gave similar results. The precision of this ranging 
technique was determined from the variation in the lidar-measured range while flying level over a flat 
surface. We chose the flight segments over the flat playa at Railroad Valley NV and calculated at the 
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standard deviation in the measured optical range for each 1 s of the 10 s. The median standard 
deviation for -140 of these 10 s segments was 1.2 m. The range accuracy depends on the timing 
accuracy of the instrument and can be biased by offsets in the receiver electronics. The lidar’s timing 
electronics have a timing accuracy of around 1 in 10 5 , corresponding to a range resolution limit of <0.1 m. 
We measured and corrected for timing offsets to within 8 nsec or ~1.2 m by using the reflection of the 
outgoing lidar pulses off the aircraft window as reference. This gives an overall accuracy of <3 m for 
the lidar’s optical range. 

The elevation of the reflecting surfaces varied widely during these flights. Figure 10 shows a 
histogram of ranging precisions, made over 10 s, for measurements made to the top of the marine 
stratus cloud layer over the Pacific Ocean, to the mountains surrounding Railroad Valley NV, and 
those to the flat Railroad Valley playa. The median values of the range precisions were 12.9, 25.2 and 
1.2 m respectively. The range variabilities >2 m were caused by surface elevation changes over the 10 s 
time (2-2.5 km distance) used for the comparison. For measurements to cloud tops, this includes any 
variability in the lidar pulse’s penetration into the cloud. 

Figure 10. Examples of range precision. Histograms of the standard deviation of range 
measurements are plotted for measurements made, and to the tops of marine stratus clouds 
(blue), to mountain tops (red), and over smooth topography (pink). These were computed 
from 10 s groups of the 1 s measurements. For this plot, only the data was screened to 
exclude measurements made in turns by selecting only data collected with aircraft tilt 
(pitch and roll) <30 mrad (40%-60% of data). Measurements made to the top of the 
marine stratus cloud deck during the Pacific Ocean flight (blue) showed a peaked 
histogram with a median value of 12.9 m. Measurements made over the mountains 
surrounding Railroad Valley NV (red) had considerable variability with a median value of 
25.2 m. This spread was caused by changes in the topographic height of the mountainous 
terrain during the 10 s analysis interval. During the same flight the median of 
measurements made over the flat playa (pink) was 1.2 m, slightly higher than the lidar’s 
range resolution limit. 



Standard Deviation in Range (m) 
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8. Airborne Measurement Results 

Figure 11 shows the time history of the Pacific Ocean flight. During this flight the lidar 
measurements were made to the top of the marine stratus cloud deck from aircraft altitudes that varied 
from 1.5 to 13 km. The bottom panel shows measurements of range to the cloud tops, the cloud top 
elevation above mean sea level, the measured DOD(pk,50) values, and the predicted values for 
DOD(pk,50). Although some limited penetration of the lidar pulse into the cloud is expected, our 
present analysis does not estimate the penetration distance. The DOD measurements closely follow the 
range changes, due to aircraft altitude and to roll angle, and the predictions. The top panel shows the 
retrieved values closely track the predicted 390 ppm column average value over the entire altitude 
range the top panel is from column averages based on the AVOCET in situ sensor during the spiral. 

Figure 11. Time history of the Pacific Ocean flight over a marine layer stratus cloud deck 
with altitudes varying from 1.5 to 13 km. The bottom panel shows measurements of range 
to the cloud tops (blue), the cloud top elevation above mean sea level (green), the 
measured DOD values (red), and the predicted DOD values (black). The top panel shows 
the retrieved concentration from the lidar measurements. The measurements plotted are 
averaged for 10 s. The measured DOD tracks the calculated DOD throughout the flight 
segment. The cloud elevation is precisely determined from the lidar range and aircraft 
altitude, pitch and roll angles. This allows for accurate measurement of the DOD as can be 
seen in the upper plot. The top panel shows the retrieved values closely track the predicted 
390 ppm colu mn average value over the entire altitude range the top panel is from colu mn 
averages based on the AVOCET in situ sensor during the spiral. The bottom plot shows the 
aircraft’s position vs. time. 



Remote Sens. 2014, 6 


458 


Figure 12 shows a time history of measurements for the flight over Railroad Valley NV, in a similar 
data format. This location is a salt flat surrounded by mountains, so the oval flight had small flat 
segments of topography surrounded by those with rapidly varying ground elevation. Since the range is 
also measured by the lidar, the measured and calculated DODs follow each other closely, even as they 
both vary rapidly during passes over the mountains. The retrieved concentration closely follows the 
predicted column average value of 390 ppm calculated from the AVOCET measurements. 


Figure 12. Measurements for the flight over Railroad Valley NV, with the data format, 
colors and averaging the same as in Figure 10. This location is a flat playa surrounded by 
mountains, so the flight track had rapidly varying ground elevation. The measurements to 
the playa are evident from the flat segments in the ground elevation plot. The highly 
variable range and DOD are due to the varying topography, aircraft banked turns over the 
mountains, and the increasing aircraft altitude in successive passes over the ground track. 
Since the range is also measured by the lidar, the measured and calculated DODs follow 
each other closely, even as they both vary rapidly during passes over the mountains. The 
retrieved concentration closely follows the predicted column average value of 390 ppm 
calculated from the AVOCET measurements. 
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Figure 13 shows the time history of lidar measurements for the flight over Four Comers NM. The 
lidar measured values of DOD(pk,50) closely follow the range variations for the 5.5 hour flight, as 
expected. The retrievals also show a nearly constant CO 2 concentration of ~390 ppm over the range of 
flight altitudes from 3-12 km, except for a 40 ppm spike at about 17:20 UTC. This was measured 
when the DC-8 directly above the Four Comers coal-fired power plant. 
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Figure 13. The retrieved values of range, DOD and CO 2 concentrations from the lidar for 
the Four Comers flight, with the data format, colors and averaging the same as in Figure 10. 
The retrieved values of DOD from measurements closely follow the range variation during 
the 5.5 h long flight, as expected. They also show a nearly constant CO 2 concentration of 
-390 ppm over the range of flight altitudes from 3 to 12 km, except for a 40 ppm spike at 
about 17:20 UTC measured directly above the Four Corners power plant. Note: The y-axis 
scale has been expanded to accommodate the spike. The inset shows an expanded view of 
flight segment from 17:10 to 17:50 UTC. This contains the spike in CO 2 that causes the 
higher DOD value beneath the arrow. 



Figure 14 shows the atmospheric conditions for the lidar measurements over Iowa. The photograph 
was taken from the aircraft window when over Iowa and shows an example of the cirrus and broken 
cumulus cloud layers that were near the spiral-down site. There are two plots of backscatter profile 
history on adjacent 30 min spans recorded from the lidar during the aircraft’s stair step ascent. These 
show the range resolved scattering from the cirms clouds between 8-10 km altitude, from broken 
cumulus cloud tops at about 2 km altitude, and increased scattering from the boundary layer below 2 km. 
The scattering from the ground is the solid red line at a few hundred meter altitude. 

Figure 15 shows the time history of the lidar measurements of range, DOD and C0 2 concentrations 
for the Iowa flight. The range changes are from aircraft altitude changes and the increases of range 
during the aircraft’s banked turns are evident. The retrieved DOD values closely follow the measured 
range. The measurements show starting at 20 h UTC when flying eastward over the Great Plains, there 
is a decrease in atmospheric C02, likely due to growing vegetation in the mid- west. The impact of the 
lower CO 2 concentration within the boundary layer for this flight is the most noticeable in the 
retrievals for lower flight altitudes, where low boundary in the boundary layer is a larger column 
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fraction, so the predicted column average is also lower. The predicted values for the column average 
based on the AVOCET measurements are shown in black in the top panel. Even with the presence of 
cirrus clouds and broken cumulus clouds the retrieved lidar values closely follow those from the 
in situ sensor. 

Figure 14. Atmospheric conditions for measurements over Iowa. (Top) Photograph taken 
from the aircraft window of the cirrus and broken cumulus cloud layers over the site. 
(Middle) Backscatter profile history recorded from the lidar from 22:08 to 22:38 UTC. 

The vertical axis is elevation above mean sea level and the top of the blue band is ~1 km 
below the aircraft altitude. Beneath that is the R 2 -corrected backscatter profile, following 
the color bar (linear scale) on the right. The scattering from the cirrus clouds between 8-10 
km altitude, from broken cumulus cloud tops at about 2 km altitude, and the boundary 
layer below 2 km are evident. The ground is the solid red line at a few hundred meter 
altitude. (Bottom) Similar plot for the next 30 min time period of 22:39 to 23:09 UTC. 
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Figure 15. The retrieved measurements of range, DOD and CO 2 concentrations for the 
Iowa flight, with the data format, colors and averaging the same as in Figure 10. The 
retrievals follow the algorithm shown in Figure 4. The retrieved DOD values closely 
follow the measured range. Starting at 20 h UTC when flying eastward over the Great 
Plains, there is a decrease in CO 2 due to growing vegetation in the mid-west. The impact of 
the lower CO 2 concentration within the boundary layer for this flight is the most noticeable 
in the retrievals for lower flight altitudes, where the predicted column average has a lower 
value. The predicted values for the column average based on the AVOCET measurements 
are shown in black in the top panel, and the retrieved values from the lidar follow 
them closely. 



Figure 16 shows a summary of the altitude resolved comparison of the column average retrievals 
from the lidar measurements with those measured by the AVOCET in situ sensor for all flights. The 
altitude resolved measurement statistics for the flights are also summarized in Table 3. In the plots the 
lidar measurement error bars are ±1 standard deviation for a 10 s average. Results are shown both for 
retrievals based on the MERRA atmospheric model and on the actual atmospheric conditions 
measured by the DC-8 during its the spiral maneuver. The retrievals based on the DC-8 measured 
atmosphere show a closer match to those from the AVOCET in situ sensor. The impact of the lower 
CO 2 concentration in the Iowa’s flight boundary layer is evident as the trend toward smaller CO 2 
column densities for measurements made from lower flight altitudes. Overall the lidar retrievals are in 
quite good agreement with the in situ measurements, particularly for those made at 8 km and higher 
altitudes, and for all made over Iowa. 
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Table 3. Summary of 2011 airborne measurement results. These use retrievals based on 
the MERRA atmospheric model as well as those on the atmosphere measured by the DC-8 
aircraft in the spiral. 


Location & Summary 
Statistics 

Mean 

Altitude 

(km) 

Min & 
Max 
Altitude 
(km) 

# of 
Meas 

Offline 

Single 

Wavel. 

Ave 

Signal 

(Counts) 

Avocet 

Ref. 

xco 2 

(ppm) 

Lidar 

Mean 

xco 2 

Based on 
MERRA 
atm 
(PPm) 

Lidar 

SD 

(ppm) 

Lidar 

Mean 

xco 2 

Based on 
MERRA 
atm. 
In-situ 
(PPm) 

Lidar 

Mean 

xco 2 

Based 

on 

DC-8 

atm 

(PPm) 

Mean 

Lidar 

(DC-8) 

In-situ 

(PPm) 

Pacific Ocean Low 
Cloud Top 

3.2 

3.1, 4.0 

102 

7,125 

390.9 

395.4 

4.72 

4.6 

389.8 

-1.0 

Excess noise limit =1.8 

6.4 

6.0, 7.0 

117 

5,375 

391.3 

392.7 

2.68 

1.4 

390.0 

-1.3 

60% conf. interval = 30 ppm 

7.5 

7.0, 8.0 

31 

4,625 

390.9 

392.9 

3.27 

2.1 

390.8 

-0.1 

# process; attempts: 559/656 

8.5 

8.0, 9.0 

41 

5,250 

390.6 

392.1 

2.41 

1.6 

390.2 

-0.4 

Process success rate: 85% 

9.5 

9.0, 

10.0 

68 

4,875 

390.5 

393.5 

2.43 

3.0 

391.8 

1.3 


12.7 

12.0, 

12.8 

102 

2,375 

390.5 

392.6 

4.05 

2.1 

391.7 

1.2 

RRV to Ground 

4.8 

4.5, 5.0 

51 

7,000 

392.7 

385.6 

4.59 

-7.1 

389.6 

-3.1 

Excess noise limit = 3.0 

6.4 

6.1, 7.0 

179 

7,750 

391.9 

388.6 

3.09 

-3.3 

391.3 

-0.5 

60% conf. interval =10 ppm 

8.0 

8.0, 8.2 

41 

7,625 

391.4 

390.0 

2.31 

-1.4 

392.0 

0.6 

# process; attempts: 514/632 

9.6 

9.0, 9.7 

55 

5250 

391.1 

389.1 

2.89 

-2.1 

390.8 

-0.4 

Process success rate: 81% 

11.2 

11.1, 

11.6 

53 

6,000 

390.9 

389.6 

2.76 

-1.3 

390.8 

-0.1 


12.8 

12.8, 

12.9 

84 

6,125 

390.9 

391.4 

2.85 

0.5 

391.4 

0.4 

Four Corners to Ground 

4.8 

4.0, 4.9 

643 

9,500 

391.8 

384.9 

4.63 

-6.9 

390.8 

-1.0 

Excess noise limit =1.8 

6.5 

6.4, 7.0 

96 

7,625 

391.5 

387.4 

3.62 

-4.1 

391.3 

-0.2 

60% conf. interval =10 ppm 

9.7 

9.4, 

10.0 

52 

7,000 

391.1 

389.3 

3.16 

-1.8 

391.7 

0.6 

# process; attempts: 846/971 











Process success rate: 87% 











Iowa to Ground 

3.2 

3.0, 3.9 

71 

16,000 

374.4 

375.9 

4.60 

1.6 

376.8 

2.5 

Excess noise limit =1.8 

4.6 

4.0, 5.0 

123 

9,000 

378.0 

377.9 

3.29 

0.0 

376.8 

-1.1 

60% conf. interval =10 ppm 

6.3 

6.0, 7.0 

43 

7,000 

380.0 

381.1 

2.62 

1.2 

381.1 

1.1 

# process; attempts: 536/704 

7.6 

7.0, 8.0 

47 

7,625 

380.8 

380.9 

1.72 

0.1 

381.9 

1.1 

Process success rate: 77% 

10.8 

10.1, 

11.0 

72 

6,375 

382.7 

381.9 

1.78 

-0.8 

383.2 

0.4 


12.5 

12.1, 

12.6 

81 

5,500 

383.8 

384.1 

2.19 

0.3 

384.1 

0.4 
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Figure 16. Comparison of column average retrievals from the airborne lidar measurements 
vs. altitude, for flights above the (a) Pacific Ocean, (b) Railroad Valley, (c) Four Comers 
NM and (d) Iowa. The plot shows lidar retrievals using the atmospheric profiles from the 
MERRA model (red), and those using the atmosphere sampled by the DC-8 during the 
spiral down (blue), and the column average CO 2 values from the AVOCET in situ sensor 
(black). The lidar measurement error bars are ±1 standard deviation for a 10 s average, and 
they apply to retrievals based on MERRA and DC-8 based atmospheres. The impact of the 
lower CO 2 concentration in the Iowa boundary layer is evident as the trend toward smaller 
CO 2 column densities for measurements made from lower flight altitudes. The lidar 
retrievals based on the actual atmospheric conditions measured by the DC-8 in the spiral 
show a close match to those measured by the in situ sensor. Overall the lidar retrievals 
based on the actual atmosphere are in quite good agreement with the in situ measurements, 
particularly for those made at 6 km and higher altitudes. The numerical values plotted are 
also summarized in Table 3. 



350 360 370 380 390 400 410 350 360 370 380 390 400 410 
C0 2 Cone, (ppm) CO ? Cone, (ppm) 


8. Summary 

We have adapted our previously demonstrated pulsed direct detection IPDA lidar for CO 2 [25] to 
fly on NASA’s DC-8 aircraft. The lidar measures the atmospheric backscatter profiles and shape of the 
1,572.33 nm absorption line using 250 mW average laser power, 30 wavelength samples per scan with 
300 scans per second. Our post-flight analysis estimates the lidar range and pulse energies at each 
wavelength every second. We then solve for the optimum CO 2 absorption line shape, and calculated 
the Differential Optical Depth at the line peak and the column average CO 2 concentrations. We 
compared these to radiative transfer calculations based on the HITRAN 2008 database, the 
atmospheric conditions, and the CO 2 concentrations sampled by in situ sensors. 

During August 2011 we participated in the ASCENDS science flights and report lidar 
measurements made on four flights over a variety of surface and cloud conditions near the US. These 
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included over a stratus cloud deck over the Pacific Ocean, to a dry lake bed surrounded by mountains 
in Nevada, to a desert area with a coal-fired power plant, and from the Rocky Mountains to Iowa, over 
cloud land with both cumulus and cirrus clouds. Most flights had 5-6 altitude steps to >12 km. 

Analyses show the retrievals of lidar range, CO 2 colu mn absorption, and CO 2 mixing ratio worked 
well when measuring over topography with rapidly changing height and reflectivity, through thin 
clouds, between cumulus clouds, and to stratus cloud tops. The retrievals shows the decrease in 
column CO 2 due to growing vegetation when flying over Iowa cropland as well as a sudden increase in 
CO 2 concentration near a coal-fired power plant. For regions where the CO 2 concentration was 
relatively constant, the measured CO 2 absorption profile (averaged for 50 sec) matched the predicted 
profile to better than 1% RMS error. For 10 s averaging, the scatter in the retrievals was typically 
2-3 ppm. The measurement precision was limited by the linear range of the PMT detector and the 
signal photon count for the wavelengths tuned onto the absorption line. This is consistent with the 
error analysis in [25], These retrievals did not use free parameters or empirical adjustments, and >70% 
of the measurements passed screening and were used in analysis. The summary in Table 3 shows the 
differences between the mean lidar retrieved values, based on the DC-8 measured atmosphere, and the 
in situ measured CO 2 column concentrations was <1.4 ppm for all four flights for measurement 
altitudes >6 km. Although the causes of these small differences are not known, we suspect a 
combination of instrument drift, instrument calibration errors, and small errors in the spectroscopic 
parameters for the CO 2 line. 
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Appendix 

Appendix A: Lidar Wavelength Calibration and Monitoring 

The lidar’s DFB laser seed laser is swept in wavelength by using a saw-tooth voltage ramp. The 
voltage from the ramp tunes the DFB seed laser -100 pm, and it was optimized to give constant sweep 
rate and a fast flyback. The power from the seed laser is divided between the transmitter and the 
wavelength calibration and monitoring setup using a 90/10 fiber splitter. The 10% fiber tap, used for 
calibration, is further divided using a 50/50 fiber splitter between a fiber coupled CO 2 absorption cell 
and the wavelength subsystem. The ramp is repeated every 3.0 ms and the sweep and the pulse 
generation are carefully synchronized. 

A heterodyne approach is used to monitor and calibrate the wavelength sweep of the DFB seed 
laser. Each calibration scan of the sweep is performed while the CO 2 sounder is operating and each 
was completed within a few minutes. The calibration subsystem used a CW DFB laser as a optical 
local oscillator and reference marker, and its wavelength is monitored using an optical wavemeter 
with an uncertainty of ±0.2 pm (±0.3 pm at 1,500 nm). The calibration subsystem mixes the 
wavelength-scanning CW DFB laser seed’s output on a detector with ~10 MHz bandwidth. 

When the seed laser’s frequency is within a few megahertz of the local oscillator a heterodyne 
signal is generated from the detector. This beat note signal is recorded as a function of voltage ramp 
scan time and averaged. The statistical center of the -20 MHz wide burst is the seed laser’s 
wavelength at a particular time in the ramp The beat-note noise burst and time with respect to the start 
of the scan is recorded using 16 bit DAQ at a 1.25 MHz rate. The local oscillator laser is then tuned to 
a new wavelength and the process is repeated. By tuning the laser to 30 positions in the scan, evenly 
distributed in time, the wavelength of the individual pulses can be determined to -20 MHz in 
frequency. The results from the flight showed during many hours of operation, the overall wavelength 
span, and line spacing remains constant except for a slow drift in the wavelength at the start of the 
sweep. This drift of <50 MHz/h was also tracked in real time by using the CO 2 reference cell 
independently of the heterodyne measurements. 
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Appendix B: Lidar Data Processing and Screening 

We used the automated procedure shown in Figure A1 to assess and screen the CO 2 column number 
density measurements from the airborne lidar. The individual returns of the one-second records were 
first sorted by range so that cloud and surface echo pulses were separated into layers. Each layer 
consisted of echo pulses from the either the ground surface or a single layer of cloud top. Cloud tops 
were identified based on similarity in range. Once sorted by range, the echo pulse train from each layer 
could be averaged to reduce shot noise. 


Figure Al. Screening procedure used to accept or reject the 10 s averaged lidar 
measurements of XCO 2 for this paper. Here CI 60 = 60% denotes the 60% confidence 
interval and XNR denotes excess noise ratio. The most common reasons for measurement 
rejection are given in the right side of the flowchart, and typical fractions of data rejected 
at each step are shown in red. The case of insufficient averaging is caused by thick clouds 
that block the lidar path to ground. For the regions where we made cross comparisons 
between lidar and Avocet measurements, >75% of 10 s data passed all the screening criteria. 
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For the averaged data, the C0 2 mixing ratio, XC02, was then determined by adjusting it to 
minimize the the difference in the measured transmission lineshape T(X) and the model lineshape, 
TcaicA X). That is by finding the value of XCO 2 that minimized the mean square weighted sum, 

x 2 (X)= X mi) [T(A) - T calc a X)l 2 /Z W(l) (B. 1) 

Here JV(X) is a weighting function that accounts for the relative sensitivity of T(A) to changes in total column 
XC02 and shot noise, and is given by 

(B.2) 


mi) = p(a) odu)N [p(a)+b]. 
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In this expression P(X) is the number of detected photons in the return pulse with wavelength X, and 
B is the number of counts from solar background and detector dark current in each pulse. 

We also calculated a confidence interval for XCO 2 . This gives the span of values for Xco 2 within 
which x 2 does not increase by more than certain factor. For this analysis we chose to limit the span of 
X 2 to within 1.15 of x 2 m in corresponding to a 60% confidence interval, or 

CI 60 =X(X 2 = 1.15x 2 mi„W>Xco 2 )-X(x 2 = 1.15 X 2 min, X<X C02 ) 

We also defined the excess noise ratio, XNR, as 

XNR = V P VM?(k) 2 ] (B.3) 

where 

XNR(X) = V[ Pa) (7U)/r calc a Xcoi) - l) 2 ] (B.4) 

Here P(X) is the number of detected photons in the return pulse with wavelength X, and T(X) is the 
transmittance after adjusting the CO 2 concentration. 

We screened the CO 2 measurements using the process shown in Figure Al. We used the size of the 
confidence interval, and a noise estimate as the primary screening criteria. For these flights the primary 
acceptance criteria were if CI 60 < 10 ppm and the XNR < 1.8. The confidence interval size criterion 
was relaxed to 30 ppm for the cloud deck in the Pacific Ocean flight due to the weak returns from the 
cloud tops. The excess noise ratio limit was relaxed to 3.0 for the Railroad Valley flight owing to 
difficulties in modeling the mountainous terrain with MERRA. 

When processing data in 10 s averaged segments, we sometimes found anomalously low 
measurements arising from one or more bad 1 s readings. Such bad readings often occurred when 
flying directly over water bodies such as the ocean, rivers, and lakes, usually due to strong specular 
reflections driving a nonlinear response from the detector. These usually had significant excess noise 
compared to the expected shot noise. 

We excluded the bad readings (<10% of the Is data), prior to averaging by processing individual 1 
s data without averaging and using the initial screening criteria shown in Figure Al. Some outliers 
remained after screening for these criteria. These were eliminated by testing against the three 
additional conditions shown in Figure Al: (a) Fewer than 3 readings for averaging in a 10 s interval. 
Such measurements don’t average sufficiently and are susceptible to anomalies, (b) Off-line pulse 
counts of less than 3,750. Such readings are sensitive to small changes in the background level, 
(c) Lidar range < 3,750 m. Such measurements had low absorption and were sensitive to small changes 
in the instrument baseline. For all flights, >75% of the 10 s measurements passed all screening criteria. 
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